Unconventional Spin Density Waves in Dipolar Fermi Gases 
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The conventional spin density wave (SDW) phase [lj, as found in antiferromagnetic metal for 
example 2], can be described as a condensate of particle-hole pairs with zero angular momentum, 
I — 0, analogous to a condensate of particle-particle pairs in conventional superconductors. While 
many unconventional superconductors with Cooper pairs of finite I have been discovered, their 
counterparts, density waves with non-zero angular momenta, have only been hypothesized in two- 
dimensional electron systems 0]. Using an unbiased functional renormalization group analysis, we 
here show that spin-triplet particle-hole condensates with 1=1 emerge generically in dipolar Fermi 
gases of atoms [2| or molecules [H, [f| on optical lattice. The order parameter of these exotic SDWs 
is a vector quantity in spin space, and, moreover, is defined on lattice bonds rather than on lattice 
sites. We determine the rich quantum phase diagram of dipolar fermions at half-filling as a function 
of the dipolar orientation, and discuss how these SDWs arise amidst competition with superfluid 
and charge density wave phases. 

PACS numbers: 



The advent of ultra-cold atomic and molecular gases 
has opened new avenues to study many-body physics. 
One of the central subjects of condensed matter physics 
is quantum magnetism, a phenomenon that has in- 
trigued scientists for decades. A quintessential example 
is the square-lattice Fermi-Hubbard model at half-filling, 
which, even at weak coupling, exhibits SDW, with the 
well-known checkerboard pattern depicted in Fig. [Ha). 
Interestingly, the theory of such antiferromagnetic or- 
der can be cast completely analogous to that of s-wave 
superconductivity - with condensation of particle-hole 
pairs corresponding to condensation of Cooper pairs in 
the BCS theory [3J. This analogy is quite robust and 
extends to particle-hole pairing with higher angular mo- 
mentum, predicting the existence of a whole array of 
SDW states, SDW„, where the label v = s,p,d,.. indi- 
cates I = 0, 1,2, .. respectively The familiar antifer- 
romagnetic order of Fig.QJa) constitutes SDW S , whereas 
SDWp is the particle-hole analogue of the spin-triplet 
p-wave superconductors and superfluid 3 He In the 
charge sector, a similar analogy predicts the existence of 
generalized i/-wave charge density wave (CDW^). 

While several candidate systems have been discussed 
for spin-singlet charge density waves (CDW^s) 0, HI, 
up to now the realization of spin-triplet SDW„^ S has 
remained elusive. These exotic states display a modula- 
tion of spin vector <S defined on the bonds of the lattice, 
such as the checkerboard pattern of red arrows depicted 
in Figs, [ljc) and [lid), as opposed to modulation of on- 
site variables in conventional SDW S and CDW S shown in 
Figs. Q^a) andQJb). Specifically, the order parameter of 
SDWp, featuring particle- hole pairs with p y -orbital sym- 
metry, is related to S v — (at aP^gCij^) for relative co- 
ordinate Yi — Vj = y (all distances are in units of the 



lattice constant throughout this paper, and repeated in- 
dices are summed over). Here a 11 with r\ £ {x,y,z} are 
the Pauli matrices, and a^ a is the fermionic annihila- 
tion (creation) operators for pseudo-spin a £ a t 
site i. The SDW s +d shown in Fig. (Ud) contains an ex- 
tended s-wave and a d xy wave component, and its or- 
der parameter is defined similarly with correspond- 
ing to diagonally opposite sites, e.g., — rj = x + y. 
In contrast, the conventional SDW S and CDW S are de- 
scribed by on-site order parameters (a] ^cr^Si,^) and 

(fit) = E a ("i,«) = EaW.cAa) respectively. 

The key insight of this paper is that fermions in a 2D 
lattice with dominant dipole-dipole interaction have the 
right ingredients to stabilize p- and d-wave SDWs. They 
emerge between the CDW and the BCS regime in the 
phase diagram, as a result of the competition between 
the short-ranged inter-atomic and the anisotropic long- 
ranged dipolar interaction. 

In a new generation of experiments, ultra-cold gases 
of dipolar fermions have become accessible in the quan- 
tum degenerate limit. Fermionic atoms of dysprosium 
161, with a large magnetic moment of 10 Bohr magne- 
ton, have been successfully trapped and cooled well below 
quantum degeneracy j4|. The fermionic polar molecule 
40 K 87 Rb has been cooled near quantum degeneracy Q 
and loaded into optical lattices. Recently, the formation 
of ultra-cold fermionic Feshbach molecules of 23 Na 40 K 
has been achieved @. On the theory side, many body 
physics of single-species (spinless) dipolar Fermi gases 
have been explored by many groups. Numerous quantum 
phases are predicted: charge density wave [9l-ll2l|. p-wave 
superfluid |13l4l7|. liquid crystalline [18l42C| . supersolid 
|23| . and bond-order solid Q. 

Here we consider a two-component (pseudo-spin 1/2) 



V d <£U 



v d ^u 




v d ^u 






t u i 








61 = 0,^ = 



v d ^u 



V r = %f(9,<f>) 



-x 



FIG. 1: Schematic spin and charge order for pseudo-spin 1/2 dipolar fermions in two dimensions. The central 
image illustrates the onsite interaction U between opposite spins (red arrows) and the dipolar interaction V T , here assumed 
to be spin-independent. The characteristic scale for V T is Vd = d 2 /a 3 , with a the lattice spacing. V T depends also sensitively 
on the orientation of the dipoles labelled by angles (6,<f>). (a) The conventional antiferromagnetic spin density wave (SDW S ). 
(b) Checkerboard charge density wave (CDW S ), where "charge" is defined as the total density, (c) An example of p-wave spin 
density waves (SDW P ) with modulation of y bond variables, (d) An example of mixed (extended) s- and d-wave spin density 
waves (SDW s+1 i). Red arrows in (c) and (d) indicate the direction of the spin vector <S defined on the bonds (yellow ellipsoids). 



dipolar Fermi gas 24 26] in an optical square lattice at 
half filling. The two pseudo-spin states can be two hy- 
perfine states of Dy atoms, or two rovibrational states 
of KRb molecules. This provides a tunable platform for 
quantum simulation of interacting fermions with long- 



2ll . [22fl ) beyond the Fermi- Hubbard 



range interactions 
model. The system is described by the Hamiltonian 



2 / j ^H,(r^i 

i.cr 



fnhj. (1) 



The lattice is aligned along the x- and y-directions, with 
nearest neighbor hopping t and on-site interaction U . U 
contains contributions from the bare short range inter- 
action, and the on-site dipolar interaction V^ , defined 
below. We assume that all dipoles are aligned in the 
same direction d = dd = (d, 9, 4>) by an external mag- 
netic (or electric) field. In general, the off-site dipole- 
dipole interaction can be decomposed into equal- and 
unequal-spin components, labeled by || and _L, respec- 



tively, 
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[1 - 3(f • d) 2 ]/r 3 . We will mostly assume V d L (d) = 

vj'(d) = Vd(d), as in Eq. [T] which arises naturally when 
the two states are associated with the same hypcrfine 
manifold. The V^{d) ^ vj{d) case will be discussed 
briefly further down. 



To give a heuristic argument about possible orders of 
the system, we consider a simplified version of model 
(HJ) retaining only the nearest and next-nearest neigh- 
bor dipolar interactions, denoted VWjV) and V x + y respec- 
tively, see Fig.[T] First, for d = z, dipolar interactions are 
purely repulsive. For U 3> Vd, the Hamiltonian reduces 
to the Fermi-Hubbard model, implying a ground state 
with SDW S order at half-filling, Fig. Ufa). For U < V d , 
the dipolar energy is reduced by placing same-spins on 
diagonally opposite sites, while opposite spins share the 
same site with only a small energy cost U. This implies a 
checkerboard modulation of the total density m = {hi}, 
i.e. CDW S , shown in Fig. (lib). 

As d is tilted away from z towards the ^-direction, 
there exists a region of tilting direction for which the 
nearest neighbor interaction V x becomes attractive while 
Vg and Vx+y remain repulsive. For instance, for <f> = 0, 
this region is bounded by two critical values of 9: $ c i = 
cos" 1 (72/3) w 35° and tf c2 = sin^y^) « 54°. In 
the simpler case of spinless dipolar fermions, a checker- 
board bond order solid is formed [§| in this region. Then 
it is plausible that for the spin 1/2 case, unconventional 
SDWs of non-s wave symmetry may be stabilized by 
interaction-induced correlated hopping either along the 
x, y, or the diagonal x + y direction. The spatial sym- 
metry of these SDWs depends on the value of 4>. This 
scenario is illustrated in Figs. [He) andQJd). 

Finally, for large dipole tilting angles, e.g., 9 > i) C 2 



FIG. 2: Phase diagram of dipolar fermions on a square lattice at half filling obtained from FRG. It is shown on 
the surface of a sphere as a function of the dipole orientation angle 9 and <j) f° r fixed interactions (a) Vd = 0.5, U = 0.1; and 
(b) Vd — 0.5, U — 0.5 in units of t. Unconventional SDWs (SDW P and SDW S+C j) are sandwiched between the CDW and BCS 
superfiuid phase. The FRG eigenwavefunction corresponding to a representative point (marked by "*") in each phase is shown 
in the k x — k y polar plots with matching colors. Note that FRG predicts a mixed s + d-wave (rather than pure d-wave) SDW. 
As the on-site interaction is increased from U = 0.1 in (a) to U = 0.5 in (b), the SDWs+d phase expands and squeezes out the 
neighboring phases. The s-wave component of SDW s +d increases with U, while the d-wave component diminishes. 



for <f> — 0, the dominant dipolar interaction is attractive. 
The leading instability is towards formation of Cooper 
pairs. Again, the precise orbital symmetry of the BCS 
phase is determined by the value of <j). 

We now determine the phase diagram of these intricate 
competing orders in the weak coupling limit, {£/, Vd} < t. 
In particular, we prove the existence of unconventional 
SDW phases for intermediate tilting angles. We use 
the functional renormalization group (FRG) technique, 
which takes an unbiased approach (without any a priori 
guess) to isolate the most dominant instability among all 

possible orders d,i3-[3Q!]. The FRG used here is an SU ( 2 ) 
symmetric version of that previously applied to treat 
spinless dipolar fermions The key ingredients of the 
FRG calculation are: (1) Derive and solve the renormal- 
ization group flows for the generalized four-point vertex 
between unequal spins, W i " L (ki, k2, IC3), where ^2^3.4) 
are incoming (outgoing) momenta in the vicinity of the 
non-interacting Fermi surface, satisfying momentum con- 
servation ki +k2 = k3 + k4, and I is the renormalization 
group flow parameter. The flow of equal spin vertex, tlf , 
is related to that of via the spin- rotation symmetry of 
H . (2) We project out the interaction channels of interest 
at each RG step, 

W, CDW (k 1 ,k a )= (2-X)^(k 1 ,k 2 ,k 1 +Q), 
Wf DW (k l5 k 2 )= -^-(ki.ka.ki + Q), 
W ; BCS (k!,k 2 ) = ^(ki.-ki.ka.-ka), 

where the exchange operator X interchanges the incom- 



ing momenta. (3) Finally we identify the most dominant 
instability of the Fermi surface from the most divergent 
eigenvalue of the interaction matrix. The correspond- 
ing eigenvector provides information about the orbital 
symmetry of the incipient order parameter. The phase 
diagram is shown in Fig. [5J 

The phase diagram displays three types of phases: 
CDW, SDW, and BCS superfiuid. We first focus on the 
case U < Vd in the vicinity of <j> = as shown in Fig. [2a). 
Consistent with our heuristic argument above, FRG con- 
firms a checkerboard CDW (CDW S ) for small 9, and a 
spin-triplet, p-wave BCS (BCS P ) superfiuid at large 9. 
For the intermediate regime, roughly between $ cl and 
•&c2 : the flow for the SDW channel diverges rapidly, domi- 
nating over the CDW and BCS instabilities on either side. 
The SDW phase shows p-wave orbital symmetry, i.e. the 
eigenvector of the SDW p phase (shown in Fig. [5]) is essen- 
tially of the form sin k y . This admits an interpretation of 
SDWp as a particle-hole analog of triplet superconduc- 
tivity/superfluidity within Nayak's classification [3] for 
generalized SDWj,. The SDW P phase found here corre- 
sponds to the class with (a^(k+ Q)a l g(k)) = S(k) -cr a jj, 
by identifying S(k) oc ssinfcj, where Q = (±7r, ±7r). The 
position space representation implies the checkerboard 
pattern of hopping amplitudes, (al a a,j^); Tj — r.; = y, 
depicted in the schematic of Fig. [TJc) . 

Additional unconventional orders with £ ^ occur in 
the vicinity of <f> = 45°, where the nearest- neighbor inter- 
action along the lattice vectors x and y is nearly equal. 
FRG predicts three more phases, CDW s+ d, SDW s+( j, 
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FIG. 3: Mean field results for the SDW P phase, (a) 

Preferred direction of the spin polarization vector S as a 
function of the ratio V} /Vj 1 ■ It is along the z-direction for 
V/ < V} and lies in the x-y plane for V/ > V} . (b) The 
energy difference between the mean field states with 5 Z and 
S x order. The inset shows the magnitude of the correspond- 
ing order parameter. The parameters are U = 0.1, V} = 1.0, 
9 = 47°, = 0. All energies are in units of t. (c) The mean 
field energy gap of the SDW P phase, in units of the Fermi 
energy, as a function of the inter-spin long-range dipolar in- 
teraction Vj~ for 9 = 47° and = 0. For the SU(2) symmetric 
case plotted in (c), the energy gaps for different vector polar- 
izations are degenerate. 



and BCS^-i-rf, all of which contain a d^-wave as well as 
s-wave components. The contributions of the isotropic 
s-wave, extended s-wave, and d-wave components are 
inferred by fitting the FRG wavefunctions using func- 
tion Co + ci cos k x cos k v + c-i sin k x sin k y , with {co , c\ , C2 } 
as fitting parameters. As a general trend, for increas- 
ing 9, the magnitude of isotropic s-wave cq reduces, 
while the magnitudes of c\ and C2 are comparable and 
increase. The CDW s +d phase can be viewed as the 
natural continuation of the CDW S as c\ and C2 be- 
come appreciable. The two representative points shown 
in Fig. GHa) for the SDW s+(i and BCS s+d are fit by 
0.05 — 0.16 cos k x cos k y — 0.18sinfca; s'mk y and 0.01 + 
0.23 cos k x cos k y — 0.19 sin k x sin k y , respectively. Since 
cq is small, the real space modulation pattern for such 



SDWg+d takes the form of Fig QJd) : atoms delocalize 
across a plaquette, in the diagonal direction perpendic- 
ular to the dipole tilting direction. In contrast to the 
triplet BCSp phase at small <fi, the BCS s+ d phase is a su- 
perfluid of singlet Cooper pairs with mixed orbital sym- 
metry, £ = 0,2. 

Next we illustrate how the phase diagram changes 
as the model approaches the repulsive Fermi-Hubbard 
model (U > Q,V d = 0). We calculate the FRG flows 
for increased on-site interaction, U = 0.5, while keep- 
ing Vd fixed at 0.5. The phase diagram is shown in 
Fig. [2fb). Since the on-site interaction U favors an- 
tiferromagnetism, the SDW P phase shrinks, while the 
SDWg+d phase extends to cover a broader region, in- 
cluding that previously occupied by BCS s +<j. Note that 
the d-wave component of SDW, even though diminished, 
is always present since the dipole interaction is kept fi- 
nite. When U is further increased such that U ^> Vd, 
only the isotropic component (cq) will survive, indicat- 
ing the SDW s+( j crosses over to SDW S , the conventional 
antiferromagnetic ordering of spins in Fig. HJa). 

To corroborate the FRG prediction of the unconven- 
tional spin density waves, we use self-consistent mean 
field theory. For a square lattice of finite size N x N, 
we impose periodic boundary conditions and retain the 
dipole interactions up to a distance of 12 lattice con- 
stants. We define the various normal and anomalous av- 
erages, Pi,a,j, a > = (tij^a^a) and TTH^ai = (<ii j<T aj j<T i). 
The corresponding mean field Hamiltonian is solved self- 
consistently by starting from an initial guess of the gen- 
eralized density matrix, and iterating until desired con- 
vergence is reached. At each step the chemical poten- 
tial is tuned to maintain half filling. The results are 
checked to be size-independent by varying N > 24. In 
Fig. [3J N is set to 28. Although mean field results are 
only suggestive, they provide an independent confirma- 
tion of the FRG results and unveil the real space pat- 
terns of <S directly in the SDW,, phases. They can also 
be used to investigate the direction of S for the general- 
ized model with (d) ^ v} (d) . We search for uncon- 
ventional SDW phases with homogeneous spin density, 
n-i.a = 1/2- In and around the SDW P region predicted 
by FRG, we indeed find solutions with order parameter 
S V = W,,a^'aB° j 3,p)i r j ~ v i = V- Further, the mean field 
energy for S x order is identical to that for S y and S z , due 

to the SU(2) symmetry of H imposed by Vj-(d) = vj (d). 
This degeneracy is lifted for ^ V d ll (d). In Fig. H 

we compare the mean-field energies of the SDW P solu- 
tion with order parameter S z and S x . The z{x)- polar- 
ized order S z {S x ) is energetically favored for vj > Vj- 

(vj < Vj-). However, the degeneracy between S x and 
S y remains. The mean field results support our inter- 
pretation of the SDWp order as schematically shown in 
Fig. [He). A similar analysis can be performed for the 
SDWg+d phase. 
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In conclusion, we have established the emergence of 
unconventional spin density wave orders, SDW P and 
SDW s+( ;, along with other exotic phases with non-zero 
angular momentum, within ultra-cold spin-1/2 dipolar 
fermions on the square lattice. These phases occupy 
a sizable region of the phase diagram mapped out via 
the functional renormalization group approach. Further- 
more, the self-consistent mean field estimation for the 
energy gaps of SDW P , shown in Fig. EJc), indicates a 
critical temperature T c ~ 0.08I> for Vd/t = 1.5. Consid- 
ering the currently reported temperature of degenerate 
dysprosium T ss 0.2Tp this suggests experimental 
accessibility of these emergent phases in the near future. 
Our study thus articulates the dipolar Fermi system as 
an intriguing and novel test bed for exotic many-body ef- 
fects. They provide a fresh perspective on systems with 
competing orders. 
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